Supersolid phase of hardcore bosons on triangular lattice. 
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We study properties of the supersolid phase observed for hardcore bosons on the triangular lattice 
near half-integer filling factor, and the phase diagram of the system at finite temperature. We find 
that the solid order is always of the (2m, — m', — m') with m changing discontinuously from positive 
to negative values at half-filling, in contrast with phases observed for Ising spins in transverse 
magnetic field. At finite temperature we find two intersecting second-order transition lines, one in 
the 3-state Potts universality class and the other of the Kosterlitz-Thouless type. 
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Since the supersolid state of matter was introduced to 
physics nearly half a century ago and its theoretical fea- 
sibility was demonstrated, U there was a long history of 
experimental attempts to find it in Nature (mostly in 
^He, see, e.g., Ref. along with numerical simulations 
and theoretical predictions for models of interacting lat- 
tice bosons. Recent years have seen a renewed interest in 
this topic. On the one hand, lattice bosons are no longer 
in the realm of idealized models and can be now studied 
in controlled experiments with ultra-cold atoms in opti- 
cal potentials On the other hand, the non-classical 
moment of inertia observed for solid '*He samples in the 
torsional oscillator experiments by Kim and Chan re- 
mains largely a mystery. 

Hardcore bosons on triangular lattice with nearest- 
neighbor repulsion y > and hopping t > Q represent 
one of the simplest (and thus most promising from the 
experimental point of view) models displaying a super- 
solid phase in an extended region of the phase diagram. 
The model Hamiltonian is given by: 
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Here 6j is the bosonic creation operator, hi = h\ bi , and 
/J, is the chemical potential. A triangular lattice of iV = 
LxL sites, with periodic boundary conditions is assumed. 
The alternative formulation of in terms of quantum 
spin- 1/2 variables Si, namely 



H = -2t 



{ij) 



(2) 

provides a useful mapping to the XXZ-magnet. The su- 
perfluid state of Eq. ^ for t » V corresponds to the 
XY-ferromagnetic state of Eq. Q , while the solid state of 
bosons is equivalent to magnetic order in the z-direction. 
At half-integer filling factor, 7i(/i = 3V) — 1/2, the model 
has an exact particle-hole symmetry. 

A robust confirmation of early mean-field predictions 
of a supersoUd phase in the ground state of (Q), 5] 
was obtained by means of Green function Monte Carlo 
(GFMC) simulations, (el The supersolid phases identified 
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FIG. 1: Schematic phase diagram of Eq. Q near half-integer 
filling factor. 



in that study for densities away from half-filling (i.e., 
for n/V > 3 and fi/V < 3), can be viewed as solids, 
with filling factors v — 2/3 and v = 1/3, doped with 
holes and particles respectively. In what follows, we 
denote them as supersolids A and B. Density correla- 
tions in A and B have \/3 x \/3 ordering with the wave 
vector Q = (47r/3, 0). In A and B the average occu- 
pation numbers on three consecutive sites along any of 
the principal axes follow the sequence (—2m, m' , m!) and 
(2m, —m', —m') respectively (it is conventional to count 
densities from 1/2 to make connection with the magne- 
tization in the spin language, m^ = — 1/2), see Fig. ^ 

The model has been investigated in a series 
of recent papers, making use of advanced numerical 
techniques.!^ 0, ^3 The proposed zero-temperature 
phase diagram is similar to that of Ref. , with the no- 
table addition of a quantum superfluid-supersolid phase 
transition at n=l/2 and t/V « 0.115 and the stable su- 
persolid state persisting for smaller values of t/V. In 
Ref. the system was thought to remain a disordered 
superfluid for arbitrary t/V. The discrepancy can be at- 
tributed to known limitations of the GFMC method. 

Based on field-theoretic, exact diagonalization, and 
other arguments, Ref. [9| hints at the possibility of the 
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FIG. 2: (Color online). Probability distributions P(n+) 
for different system sizes and temperatures at ^/V — 3 and 
t/V = 0.1. 



FIG. 3; (Color online). Probability distributions P(n+) 
for different system sizes and temperatures at fi/V — 3 and 
t/V = 0.05. 



(m,0, —m) density order in the ground state at n = 1/2 
(state C). These considerations involved, in particular, an 
analogy between the properties of Eq. (0), and those of 
the Ising antiferromagnet on the triangular lattice, in the 
presence of a transverse magnetic field [ll|. If true, there 
should exist quantum A — C and C — B phase transitions 
away from half-filling and three finite-temperature tran- 
sitions of the Kosterlitz-Thouless (KT) type. Though 
Ref. 8] finds that the ground state is of the A oi B type, 
it makes similar predictions for the finite temperature 
phase diagram at n = 1/2 which follow from the assump- 
tion that spontaneous symmetry breaking between A, B, 
and their lattice translations is described by the six-clock 
model 

In what follows, we provide strong evidence that the 
supersolid state at half-filling is always of either the A 
or B type. Our data suggest that there is a discontinu- 
ous transition from A to B at fj, — 3V similar to the I- 
order phase transition (driven by the large energy of the 
A — B domain walls). What makes it special is the exact 
particle-hole symmetry; structure factor, superfluid den- 
sity, and energy remain continuous functions of fi through 
the transition line. For the supersolid A (or B) with the 
three-fold degenerate ground state, one expects to see the 
normal-superfiuid KT and the solid-liquid 3-state Potts 
transitions, as temperature is increased. Moreover, the 
KT and Potts transitions are independent of each other 
and for n ^ 1/2 intersect on the phase diagram. The 
failure of the mean-field description and analogies with 
the transverse-field Ising model to predict the supersolid 
structure at n = 1/2 can be traced back to the U(l)- 
symmetry of Eqs. ^ and as noticed in For 
example, the (1, 0, —1) state can not be the true ground- 
state at finite t in the limit oi t/V ^ simply because 
it does not respect the particle conservation law. 

We use the worm-algorithm Monte Carlo scheme in 
the lattice path-integral representation to simulate 



Eq. 1^. Since the structure factor 
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does not distinguish between supersolids A, B, C, we 
adopt the following strategy: for each system configu- 
ration, we compute the distribution of time-averaged oc- 
cupation numbers, fik = l3~^ dr hk{T), and use it to 
determine the fraction of sites with fik > 1/2 



N 



N- 



-i^0(nfc-l/2), 



(4) 



fe=i 



where 9{x) is the Heviside function. A, B, C density 
structures correspond to n\ = 2/3, = 0, and = 
1/3. Finite systems are characterized by broad proba- 
bility distributions P{ri^)^ and the formation of differ- 
ent solid orders can be seen as the development of sharp 
peaks, as the thermodynamic limit is approached. 

In Fig.[21we show the evolution of the P(n+) distribu- 
tion for the half- filled system at V/t = 10, i.e., close to the 
superfluid-supersolid transition point, estimated [1. lol ITol| 
at V/t « 8.5. The distribution is peaked at n+ = in 
the smallest system considered (i=6), but, as the sys- 
tem size is increased, the weight is shifted toward the 
wings of the distribution. For L=18, there are already 
three peaks with comparable height. Finally, in the L=24 
system we observe only two peaks corresponding to the 
supersolid phases A and B. Though the probability den- 
sity between the peaks is still measurable, the dynamics 
of the algorithm becomes very slow; it typically takes 
millions of Monte Carlo sweeps, in order for the system 
to make a transition from the A to the B structure and 
vice versa. We have explicitly verified that configura- 
tions with « 2/3 and n+ w 1/3 have density or- 
ders depicted as in Fig. with a large contrast in den- 
sity between sublattices. We have also checked that the 
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FIG. 4: (Color online). Superfluid density in the vicinity 
of the KT transition for t/V = 0.1 and fi/V = 2.74. The 
solid line is the thermodynamic curve calculated using Eq. @ 
with K(r) deduced from the plot shown in the inset. Inset: 
solutions of the Eq. ^ for different pairs of system sizes: 
1/2 = 24, Li = 12— filled circles, L2 = 48, Li = 12— open 
circles L2 = 48, Li = 24 — filled squares. The dashed line is 
the linear fit k = 1 + 1.03(rc - T)/t with Tc/t = 0.50. 
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FIG. 5: (Color online). Structure factor in the vicinity of the 
3-state Potts transition for t/V = 0.1 and fi/V = 2.74. Inset: 
data collapse using exact critical exponents for the 3-state 
Potts model il7| and 5 = (T - Tc)/t with Tc/t = 1.035. 



V/t — 10, L~A8, T — system spontaneously devel- 
ops either A ov B order, starting from the initial config- 
uration corresponding to the superfluid phase at V/t—b. 

In Fig. 121 we show what happens at larger values of 
V/t. Now, the central peak is already absent in relatively 
small L=\2 and L=18 systems. We thus conclude that 
the nature of the supersolid state at half-integer filling 
factor is determined by the A and B structures, for all 
values of t/V for which a supersolid phase exists. 

If spontaneous symmetry breaking of the ground state 
degeneracies is described by the six-clock model : one 
should observe three finite-temperature transitions for 



systems near half-filling, and a solid phase with algebraic 
correlations "sandwiched" between the solid and normal 
liquid phases. This prediction was made in Ref. Q for 
11 = 1/2. Since the ground state was found here to be 
only of the A or B type, and we do not see why domain 
wall energies between translated A states are the same 
as between A and B states (in fact, the Landau theory 
prediction 8, 9, 12j is that A and B states phase sepa- 
rate and have different average densities even at = BF), 
the finite temperature phase diagram should instead fea- 
ture the normal-superfiuid KT and the liquid-solid 3- 
statc Potts (for n ^ 1/2) transitions breaking U(l) and 
translation symmetry respectively. At n = 1/2 we expect 
only one liquid-solid transition. An interesting question 
is whether transition lines simply intersect, or there are 
bicritical and tricritical points and I-order lines as ob- 
served for the similar model on the square lattice |l5j |. 
We performed simulations for two representative cases, 
one for constant chemical potential fi/V ~ 2.74 (or den- 
sity n m 0.44), and the other for constant t/V — 0.1. 

In Fig.2|we show typical data for the KT transition be- 
tween the solid and supersolid phases. The transition is 
smeared by logarithmic finite-size effects, but the critical 
temperature can be still determined with good accuracy 
by utilizing the well known renormalization flow and the 
universal jump of the superfluid density, ps, at Tc- The 
data analysis is as follows :|16| we define R = iTps/2mT 
(where m — l/it is the effective mass for the triangular 
lattice) and study the finite-size scaling of the data using 
KT renormalization group equations in the integral form 



41n(L2/Li 



dt 



R2 
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The microscopic (system size independent) parameter k 
is an analytic function of temperature, and the critical 
point corresponds to i? = 1 at k = 1. For T < Tc, the 
thermodynamic curve is defined by the equation 



1/R + 111 R = k{T) , 



(6) 



with K = 1 -t- k'{Tc — T). We use different pairs of system 
sizes in Eq. (jSJ to determine the k{T) curve, and obtain 
the location of the critical point from k{Tc) = 1. The re- 
sults are shown in the inset of Fig. 0] Data collapse and 
smooth analytic behavior of k{T) proves that the transi- 
tion is indeed of the KT type. We used the same protocol 
and system sizes to determine other critical points. 

In Fig. Is] we present our data for the transition into the 
state with the long-range density order. For the three- 
fold degenerate B structure this transition is expected 
to be in the 3-state Potts universality class. The crit- 
ical exponents are known exactly |l7l |: v — 5/6, and 
f3 = 1/9. We thus perform the data collapse using 
L^f'SQ = fiSL^'") where 5 = {T - Tc)/t and Tc is the 
only fitting parameter. The result is shown in the inset of 
Fig. |5| This confirms the above-mentioned expectation, 
and establishes that there is only one transition to the 
solid phase (there are no visible finite-size effects below 
Tc). 
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FIG. 6: (Color online). Finite temperature phase diagram for 
two representative cuts: the left panel is for fixed fi/V = 2.74; 
the right panel is for fixed t/V = 0.1. The solid line indicates 
a degenerate I-order transition line between supersolids B and 
A. 



Finally, we compute the phase diagram in the 
{T/t, t/V) (at constant n/V = 2.74) and {T/t,fi/V) (at 
constant t/V — 0.1) planes and observe that KT and 
Potts transition lines form a simple cross for n 7^ 1 /2, i.e., 
the corresponding order parameter fields are not strongly 
interacting, see Fig.El The transition temperature to the 
superfluid and supersolid states in this part of the phase 
diagram is determined by the hopping amplitude. Within 
the statistical uncertainties of our calculation, KT and 
Potts transition temperatures cannot be distinguished at 



^i/v=3. 

We did not see evidence for the algebraic solid state at 
/i = 3V. The finite-size scaling for the supersolid-solid 
transition at /i = 3V is consistent with the 3-state Potts 
universality, though the data collapse is not as impressive 
as in Fig. ^ (the other alternative is the KT transition). 

It is instructive to understand why the (to, 0, —to) 
phase for the Hamiltonian |^ is not an obvious ground- 
state. At the mean-field level, C has a better ene rgy 
than A or B. For the transverse-field Ising model [ij 
the (1,0,-1) spin arrangement is obtained by orienting 
the middle spin along the magnetic field direction, i.e. 
putting it in the equal-amplitude superposition of up- 
and down-states. In bosonic language, it corresponds 
to the superposition of states with one or zero particles 
on a given site. Such a state can not be reconciled with 
the Hamiltonian which conserves the particle number. 
Any non-integer average occupation number necessarily 
involves hopping transitions to the nearest neighbor sites. 
In the (1,0,-1) structure the middle site is completely 
surrounded by the fully occupied or empty sites and thus 
can not be the ground state of the system even in the 
limit of t/V — !■ 00. The problem appears to be inherently 
quantum with no obvious solution at the mean-field level. 
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